Dynamical evolution of quasi-circular binary black hole data 
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We study the fully nonlinear dynamical evolution of binary black hole data, whose orbital pa- 
rameters are specified via the effective potential method for determining quasi-circular orbits. The 
cases studied range from the Cook-Baumgarte innermost stable circular orbit (ISCO) to signifi- 
cantly beyond that separation. In all cases we find the black holes to coalesce (as determined by 
the appearance of a common apparent horizon) in less than half an orbital period. The results of 
the numerical simulations indicate that the initial holes are not actually in quasi-circular orbits, but 
that they are in fact nearly plunging together. The dynamics of the final horizon are studied to de- 
termine physical parameters of the final black hole, such as its spin, mass, and oscillation frequency, 
revealing information about the inspiral process. We show that considerable resolution is required 
to extract accurate physical information from the final black hole formed in the merger process, and 
that the quasi-normal modes of the final hole are strongly excited in the merger process. For the 
ISCO case, by comparing physical measurements of the final black hole formed to the initial data, 
we estimate that less than 3% of the total energy is radiated in the merger process. 
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I. INTRODUCTION 

One of the most pressing problems in gravitational 
physics is the solution of the two body problem for an 
inspiralling pair of black holes (BHs) from quasi-circular 
orbit to merger. This process is expected to be a likely 
candidate for early detection by the current generation 
of gravitational wave detectors. Although considerable 
success has been made in long term, 3D evolutions of 
single distorted holes or head-on collisions ^ Ql , and hy- 
brid numerical-perturbative evolutions 0, Q , long-term 
evolutions of binary BH systems with angular momen- 
tum have seen only limited successes |E IE S 111 j up to 
timescales of a single orbit With limited evolution 
times, emphasis has been on starting from solutions to 
the constraint equations which approximate the very last 
phase of the inspiral process, close to the final plunge. 

Such a family of BH binary parameters was proposed 
by Cook ;10,j, and is based on finding minima of an ef- 
fective potential for boosted binary BHs (with throats). 
The effective potential approach was also applied later 
in for binaries constructed from puncture data, where 
similar results were found for the innermost stable circu- 
lar orbit, or "ISCO". Furthermore, recent work 
finds a very similar sequence of orbital parameters for bi- 
nary puncture BHs using an entirely different construc- 



tion based on the assumption of a "helical Kil ling vec- 
tor" (HKV). On the other hand, Refs. [Illlllll, the 
first to develop the HKV approach, solve the thin sand- 
wich equations for binary BHs (with throats), obtaining 
markedly different results for the estimate of the ISCO. 
These results agree with increasingly higher order post- 
Newtonian calculations 0| . At somewhat larger separa- 
tions, the various constructions show better agreement, 
as would be expected. 

It is critical to understand the differences between 
these data sets, in particular to investigate whether any 
of them represent late stages of physically realistic orbital 
BH mergers. The only method to study the fully non- 
linear orbital dynamics of these data sets, which are well 
into the strong-field region, is full numerical simulation. 
Recent advances, including the use of the BSSN evolu- 
tion system [TM Il9l 12(1 [2lll . new gauge conditions 1], and 
excision techniques |22l |23 , have brought about signifi- 
cant improvements in the ability of numerical relativity 
to handle binary systems. Further, by using co-rotating 
coordinates (implemented through an adjustment to the 
shift condition j2J|) we find that the BHs can effectively 
be held in place, greatly reducing dynamics in grid func- 
tions and removing complications due to the need of mov- 
ing excision regions. These improvements have allowed 
us to carry out the first systematic and fully nonlinear 
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numerical studies of the orbital dynamics of such data 
sets, constructed to be in the "near-ISCO-regime" . The 
dynamics of the entire inspiral and merger process can be 
followed to well over t « lOOM in many cases. (Distances 
and times are given in terms of the ADM mass, which is 
approximately unity for the initial data models studied 
here.) Similar techniques were recently used to evolve 
a pair of separated apparent horizons for the expected 
period of a close orbit Q • 

In this paper we apply these full nonlinear capabil- 
ities to examine the orbital dynamics of a near-ISCO 
sequence. The sequence, based on the effective poten- 
tial work of Refs. [Tol is constructed to represent 
a family of circular orbits of non-rotating, equal mass 
puncture BHs. These parameters were previously stud- 
ied 0,115 with a hybrid numerical-perturbative approach 
(see, e.g., 0jllal23)- Increasingly sophisticated over the 
years, this hybrid approach is aimed at extracting wave- 
forms if a linear regime can be reached, while only fully 
relativistic approaches are capable of studying the dy- 
namics of the complete inspiral process. The near-ISCO 
regime studied here is for initial proper spatial separa- 
tions ranging from 4.99Af (the Cook-Baumgarte ISCO) 
to 7.84M. The data are evolved at least until a common 
horizon forms, indicating that the pair has merged. In 
all cases we find that the merger occurs in roughly half 
(or less) of the expected period, based on the determina- 
tion of the angular velocity in the initial data calculation. 
These and other measures discussed below indicate that 
rather than orbiting BHs, these data sets represent nearly 
plunging BHs. 

Using techniques originally developed in pSJ g ^lsollBll 
Is^,!!^ and applied to grazing BH collisions in [MItIIsIIS^. 

we study the physical properties, such as mass and spin, 
and the dynamics of the final BH horizons formed in the 
merger process. In particular, we compare a number of 
independent measures of horizon dynamics, including the 
geometry of apparent horizons and event horizons, and 
the isolated horizon formalism 35, 3&, i37.J . which show 
consistent results for masses and spins of the final black 
holes. Comparing these results with the parameters of 
the initial data, we are able to infer some properties of 
the orbital dynamics. We find that at low resolution, the 
angular momentum of the final BH is significantly lower 
than the angular momentum of the spacetime. However, 
for the case studied in most detail at higher resolution, 
the angular momentum of the final BH formed in the 
merger is within about 20% of the angular momentum of 
the spacetime. Hence, with current resolutions, we are 
not yet able to use this to measure the angular momen- 
tum radiated. On the other hand, our highest resolu- 
tion calculations of the final BH mass indicate that less 
than 3% of Mabm is radiated, showing very good accu- 
racy of the overall energy accounting in the simulations, 
but the current resolutions are still not sufficient to use 
this measure to predict precisely the radiated energy. By 
studying the dynamics of the horizon geometry, we are 
further able to determine the oscillation frequency of the 



final BH, showing that its quasi-normal modes (QNMs) 
arc strongly excited in the merger process. 



II. INITIAL DATA 

The initial data are calculated via the puncture 
method . This variant of Brill-Lindquist initial data 
performs a conformal compactification at the interior 
asymptotically flat region of the BHs. The extrinsic cur- 
vature is given by the Bowen-York solution to the mo- 
mentum constraint, and the Hamiltonian constraint is 
solved for the conformal factor. The data are confor- 
mally flat and maximally sliced. The initial lapse is cho- 
sen to be a = 1, while the initial shift is determined by 
a quasi- rigid rotation (see below). 

Each member of the initial data sequence we consider 
has two BHs of equal mass and equal opposite momen- 
tum, without spin. Initial parameters for such systems 
in quasi-circular orbit are tabulated in |25l | , based on the 
effective potential calculations of and for punc- 
tures. In Table ^ we list the first five elements in this 
sequence. The effective potential construction singles 
out QC-0 as representing the presumed innermost sta- 
ble circular orbit, or ISCO, while other members of the 
sequence of quasi-circular orbits have initial radii progres- 
sively further out. We evolve all members of this near- 
ISCO regime to coalescence and beyond, using full nu- 
merical relativity. The measure L/M between the holes 
is the proper distance computed along a coordinate line 
from the location of one apparent horizon (AH) to the 
other. For convenience, we also provide dynamically de- 
termined coalescence times in the Table, as measured by 
the appearance of the horizons, both as computed by co- 
ordinate time and by the proper time measured by an 
observer at the origin between the two BHs. Further 
discussion of the evolution results are the subject of the 
following sections, but one can see immediately that by 
these measures, coalescence occurs well before the orbital 
time scale for all cases studied. 



III. EVOLUTION METHODS 

Formulation and gauge. The evolution of the system 
is carried out usin g t he techniques and evolution sys- 
tem ("BSSN") HI 13123 with the implementation de- 
scribed explicitly in U|. We use a hyperbolic version of 
the "l-|-log" lapse and F-driver shifts. In particular, the 
lapse a is evolved using 

dta = ~~2aK, (3.1) 

where K is the trace of the extrinsic curvature. The shift 
is implemented using 

dt(3' = FB' dtB' = dtV - 7]B' (3.2) 
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TABLE I: Elements of the initial data sequence of [2ql . list- 
ing the initial proper separations, angular momenta, and an- 
gular velocity (L/M, J/M, Mil, and corresponding orbital 
timescale Torbit). The coordinate time Iah and the central- 
observer proper time t^h* to common apparent horizon for- 
mation are given for the intermediate resolution dx = 0.06M 
runs with shift parameters (j],p,n) = (4, 1,4), corresponding 
to Figs. |5| and 13 Also shown is the proper time teij* until a 
common event horizon is found for a central observer. Due to 
computational limitations, we did not attempt to locate EHs 
for QC-2 - QC-4. 



with r] a constant dissipative parameter and F a function 
of space and time given by 



F 



3 aP 



(3.3) 



where V'bl is the time independent conformal factor for 
Brill-Lindquist data. This choice of F ensures a shift 
which is nearly static at late times, as well as ensuring 
the correct fall-off near the puncture. The factor 3/4 
guarantees that the asymptotic gauge speed associated 
with the longitudinal shift components is equal to the 
speed of light. These features were found to be important 
for long-term stable and accurate evolutions of head-on 
collisions without excision l|. Typically, we use values 
of p = 1 or 2 and n = 2 or 4. The parameter ry can be used 
to tune the rate of horizon expansion over the course of 
the evolution, with larger values leading to faster horizon 
growth. We typically use values rj e [2, 5]. 

As detailed in [2^1, we apply a quasi-rigid rotation to 
the initial values of the shift vector, 

f3%=o = ^y^^ /?^|,^o = x-^a (3.4) 

The factor of 1 /?Abl modifying the rigid rotation angular 
velocity Cl has the effect of ensuring that the shift is zero 
at the location of the punctures (a condition that is nec- 
essary if one wishes that the punctures remain at a fixed 
location in coordinate space). We find that with appro- 
priate values of in the initial values of the shift, the 
location of the horizons can be made to remain almost 
stationary on the grid. Since the size and location of 
the excision region is determined by the apparent hori- 
zon location, this allows the evolution to proceed until 
merger without requiring the excision region to move on 
the grid. This type of co-rotating shift, with dynamic 
control to adjust the co-rotating frame as the binaries 
orbit, is described in To initialize the co- rotation Cl, 



experiments were carried out to empirically determine 
values for which the horizon tended to move least. 

It is worth noting that the optimal values tended to 
be 50 — 60% of the values of specified in Table HJ As 
the latter are asymptotically determined, and further as 
slicing and shift conditions in the interior have strong ad- 
ditional influences on the evolution of the BH orbit (for 
example, the F-driver shift condition used in the evolu- 
tions has itself a co- rotating component H^D, it is diffi- 
cult to draw firm conclusions from this result. However, 
the result may be consistent with one of the main results 
of this paper, namely that the near-ISCO regime orbital 
sequences studied here coalesce rapidly, well within an 
orbital time period, as discussed in more detail below. 

Apparent horizons are determined at each time-step 
using the horizon finder AHFinderDirect described 
in 39] . As a consistency check, apparent horizons are 
also determined each 0.48Af using the earlier, accurate 
but much less efficient horizon finder AHFinder, an im- 
plementation of the method described in lift l4l] |. For 
QC-0 and QC-1, an event horizon finder |42| was used 
for further analysis, as a post-processing step. 

Numerical methods. Finite differencing is performed to 
second order accuracy, and the time stepping is carried 
out via a 3-step iterative Crank-Nicholson scheme. The 
boundaries are placed at the same coordinate location 
for all simulations reported here. The proper distance to 
the outer boundary is extended by applying a "fisheye" 
coordinate transformation to the initial data 0,0. At 
the outer boundaries, an outgoing-radiation (Sommer- 
feld) boundary condition is applied to each of the evo- 
lution variables. The boundary condition introduces a 
constraint violation, visible as a region where the con- 
straints do not converge to zero, which propagates in- 
wards from the outer boundary. To minimize the effect 
of this error on the interior dynamics, we have placed 
the outer boundaries of our grid at a distance which is 
causally disconnected (with respect to the physical speed 
of light) from the common apparent horizon at the time 
of merger. Although this does not guarantee complete 
numerical causal disconnection, it significantly reduces 
potential outer boundary effects in the neighborhood of 
the horizons during the pre-merger phase, and allows us 
to achieve second order convergence (for instance in the 
Hamiltonian constraint violation) in the region near the 
final merged horizon (Fig.^. Through convergence stud- 
ies and experiments with boundaries further away, we 
conclude that with our computational grid parameters, 
the outer boundary does not affect the results presented 
here. 

Excision. For the long-term simulations carried out in 
this paper, we excise the singularity using an extension of 
the "simple excision" techniques which have proven suc- 
cessful in evolving single BH spacetimes . We 
use the same condition on each of the fields. Rather than 
excise a cube we excise within a topologically spherical 
surface (a "lego" sphere) determined by the shape of the 
AH. A number of buffer points (at least 5) is maintained 
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FIG. 1: Convergence in Hamiltonian constraint for the QC- 
case shown at the time of formation of the common AH, 
t = 17.76M. The resolutions are dx = 0.08M, 0.06M and 
0.04M, with the evaluated constraint scaled appropriately so 
that the lines would lie on top of each other in the case of 
perfect second order convergence. The lower scale shows the 
physical distance along the x axis, while the upper scale shows 
the grid coordinate rgrid defined by 14.111 . The vertical dot- 
ted line shows the apparent horizon position. Note that al- 
though we do not obtain second order convergence near the 
outer boundary, there is good second order convergence in the 
neighborhood of the horizons. 

between the AH and the excision boundary in order to 
minimize the propagation of errors introduced by the ex- 
cision boundary condition out of the horizon. As shown 
in [23, for the gauge conditions used here this buffer 
zone is expected to be adequate to protect the BH exte- 
rior from errors coming from the excision boundary. The 
excision boundary is allowed to grow as the AH grows, 
but not to move or shrink in coordinate space ~ that is, 
points which are excised are never returned to the active 
grid. We find that our runs are second order convergent 
even in areas neighboring this excision boundary. Even- 
tual loss of convergence at late times, after merger, is 
entirely due to outer boundary effects. 



For each model, runs were performed with uniform grid 
resolutions of dx — O.OSAf and dx = 0.06M on grids of 
384 X 384 X 192 and 512 x 512 x 256 points respectively. 
For 3-level convergence tests, certain models were also 
run on 640 x 640 x 320 [dx = 0.048M) and 768 x 768 x 384 
{dx = 0.040A/), though especially the latter are at the 
limit of what even large-scale clusters can sustain. 

To allow the outer boundary to be farther from the 
strong-field region, the (Cartesian) grid was placed non- 
uniformly in the physical coordinates: the physical radius 
r of each grid point was related to the grid radius Tgiid 
by a "transition fish-eye transformation" (Q) 

X In ^cosh Il£lfLJlZ!2.^ _ In ^cosh 1^12^ — !-2.^ 

(4.1) 

where the parameters tq = 5.5M and s = 0.8M are the 
radius and width of the transition region, and a = 4 gives 
the physical scale factor outside the transition region. 
This placed the grid boundaries rgiid = 15.48M at a 
physical distance r = 45.42M from the origin. 

Fig- n displays a cross section of the Hamiltonian con- 
straint for QC-0 dX t — 17.76Af for the three resolu- 
tions dx = 0.08M,0.06M,0.04M, scaled to lie on top 
of each other if the code is converging at second order. 
At this time the common horizon has just appeared at 
the highest resolution. As is typical for these runs for 
which the outer boundary condition is not constraint pre- 
serving, a non-convergent effect marches inward from the 
outer boundary. The region disconnected from the outer 
boundary, however, shows second order convergence, in- 
cluding points in the immediate neighborhood of the ex- 
cision region. Significantly, convergence is maintained in 
the region of common horizon formation at the time at 
which it first forms. 



V. RESULTS 



IV. SIMULATION DOMAIN 

For each dataset a number of test runs on a small grid 
and with low resolution were performed in order to de- 
termine the approximate time of common horizon forma- 
tion. These times were used to construct grids with outer 
boundaries causally disconnected from the common hori- 
zon at merger time (distance to the outer boundary mea- 
sured by the physical distance '63|). Initial datasets were 
set up with punctures along the y-axis and momenta in 
the a;-direction. The z-reflection symmetry allowed the 
evolutions to be carried out on the positive z portion of 
the grid, with a symmetry along the z = plane imple- 
mented through the use of ghost-zones whose data where 
reflected from interior grid points. 



Evolutions to merger were carried out for each of the 
QC-0 to QC-4 initial data models, at the 384 and 512 
grid sizes described above (a complete set of 640 and 
768 runs was not possible due to computational resource 
limitations). We study the orbital dynamics and merger 
times for all five members of this sequence. Evolutions 
were typically carried out to just beyond the formation of 
a common apparent horizon (AH), though it is possible 
with careful parameter choices to carry out the evolu- 
tions for significantly longer, often to beyond t = lOOA/. 
At late times, however, contamination from the outer 
boundaries can become a significant source of error. 

Merger Times. The actual time at which the com- 
mon AH forms for each model is plotted in Fig.|21for the 
resolution dx = 0.06M. The coordinate time at which 
a common AH forms is clearly not a physical invariant. 
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FIG. 2; The time to appearance of a common AH for each of 
the QC models. Filled circles indicate the results of numer- 
ical experiments using dx — 0.06, from initial proper spatial 
separations indicated along the a:-axis. The upper line indi- 
cates the expected orbital period, based on the initial angular 
velocity (Table and assuming a Newtonian circular orbit. 
Empty circles indicate the fraction of an orbit before common 
AH formation. The error bars show ± the difference between 
dx — 0.06 and dx — 0.08 evolutions; the effects of variations in 
gauge parameters and outer-boundary position are generally 
smaller than these. 

and depends not only on gauge effects, but also on com- 
putational parameters. For each point in the near ISCO 
QC sequence we have therefore carried out a large num- 
ber of simulations, varying not only grid resolution and 
boundary location, but also gauge conditions, within the 
class of conditions discussed above. We find that the re- 
sults are only weakly dependent on these effects; based 
on these tests we have provided approximate error bars 
in Fig. 121 that span the general range of results we ob- 
tained over all parameter variations. Furthermore, we 
experimented with 4th order spatial and 3rd order tem- 
poral differencing, which improves the accuracy of the 
lower resolution results, but which does not change the 
essential results presented here. The variations in timing, 
over all gauges and computational parameters tested, are 
typically less than 10%, a result that was typical of ear- 
lier studies 0, where, for example, rather different 
lapse conditions were also tested with minimal changes 
in coalescence time results. 

For reference, the expected orbital periods based on the 
initial data estimate of the angular velocity from Table HI 
are plotted as square points (the top- most curve). For 
each of the models studied, a common AH appeared in 
less than half of the expected time for an orbit. This 
is true even for the further separated cases, well beyond 
the effective potential ISCO, where dynamics would be 
expected to better approximate a quasi-circular orbit. 

To gain insight into the physical interpretation of the 
merger times, the values obtained from full 3D evolu- 
tions of this near-ISCO sequence can be compared to the 
merger times for head-on collisions, with each puncture 
starting with zero linear momentum. In Ref. |44| it was 
shown that there is a remarkable agreement between the 



Newtonian free fall time, and the fully relativistic com- 
putation of the time it takes for a common AH to form 
for BHs colliding head-on. To do this comparison we car- 
ried out a set of 3D simulations, evolving initial data for 
head-on collisions of BHs represented by punctures ini- 
tially at the same distances of the QC sequence studied 
in Tabled but without boosts. We used the same lapse 
condition as described above for evolving the orbital con- 
figurations, while the shift condition was slightly modi- 
fied lU- The gauge parameter choices were p = 1, n = 4 
and rj € [2.8, 3.5], with the larger values of rj being used 
for the larger separations. We do not expect the small 
difference in shift condition to affect the results signifi- 
cantly. In Fig. 131 the results of the QC simulations are 
compared with both Newtonian infall times (dashed line) 
and relativistic head-on collisions (filled diamonds) . The 
latter two results lie close together on essentially parallel 
curves, agreeing with the results of |44|. and providing a 
baseline for coalescence times of purely plunging, head-on 
collisions. 

The lower curve in Fig. (21 indicates the proper time to 
common AH formation, tah central, as measured by an 
observer at the origin of the coordinate system. As this 
is a symmetry point of the problem, this proper time is 
invariantly defined and easily calculated by integrating 
the lapse, and removes most of the potential coordinate 
ambiguity that might be present in the other time deter- 
minations. This proper time tah central is lower than the 
coordinate time because it is measured in a region where 
the lapse is somewhat collapsed. 

It is important to point out that the results obtained 
for the full near-ISCO QC sequence are basically parallel 
to the curves for head-on collisions. Further, the proper 
time measured by observers at the origin, shown by the 
squares, are also consistent with these results, strength- 
ening the argument that gauge effects are not responsi- 
ble for the trends. Taken together, these results suggest 
that for the near ISCO sequence studied here, the BH's 
are not in fact in quasi-circular orbits, but are rather in 
an extended plunge towards coalescence. Even at nearly 
twice the distance of the ISCO (QC-4), there is no clear 
trend towards an orbit. On the contrary, as shown in 
Fig. 12 the configurations merge in progressively smaller 
fractions of an orbit. 

Studies of Horizon Physics. Much can be learned 
about the physics of the merger process, and of the fi- 
nal BH formed, through a detailed study of the horizons 
(see, e.g. [1 B IH Hfl HI)- We begin by studying the 
masses of the AHs as a function of time for each member 
of the sequence. The extent to which horizon masses are 
consistent with other quantities in the spacetime provides 
a meaningful and physical measure of the overall accu- 
racy of the simulation. As a measure of the BH mass, we 
study the irreducible mass Mi„ = ^Area/167r. In Fig. ^ 
we plot Mirr for each model in the near-ISCO sequence 
over the course of evolution. Initial individual horizon 
masses of approximately 0.5Madm are found to remain 
essentially constant, until the formation of a common 
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FIG. 3: The time to common AH formation is compared with 
various head-on coUision results. The upper hne (filled circles) 
repeats the curve of Fig. |21 showing common AH formation 
times for each of the QC models. The dashed line indicates 
the Newtonian collision time for a pair of particles falling to- 
gether from the given distance with zero initial velocity. The 
corresponding relativistic head-on collisions were simulated 
and give the results plotted as filled diamonds. Finally, the 
lower curve (open boxes) refers to the same QC simulations 
as the upper (filled circle) curve, but uses the proper time at 
the origin at the time of common AH appearance (see text). 



horizon. After merger, we find that in aU cases the re- 
sults are consistent with the ADM mass, and again hold 
a rather constant value in the subsequent evolution. Er- 
rors in the calculated mass are greater for the late merg- 
ing cases, but even in the worst cases the masses do not 
exceed the ADM mass by more than about 5%. Drift in 
the final horizon mass is quite small, even at late times 
when errors accumulate and outer boundaries might be 
expected to have an influence. Results from the two dif- 
ferent apparent-horizon finders are very similar, provid- 
ing a further consistency check. 

Figure. E^b) shows that these results are only weakly 
resolution-dependent. At higher resolutions, we see 
slightly smaller horizon masses, somewhat later merger 
times, and a much smaller upward drift of the common- 
horizon mass. Unfortunately, these results do not show 
quantitative convergence. This appears to be due to in- 
sufficient resolution, particularly in our coarsest (dx — 
O.OSAf) models. (Comparing this result to the good 
second order convergence of the Hamiltonian constraint 
at comparable resolutions (Fig. it would appear 
that the apparent-horizon diagnostics require higher res- 
olution than the Hamiltonian constraint to reach the 
quantitative-convergence regime.) 

Fig-EIb) shows that the individual-horizon masses are 
slightly greater than |Madm- This is not surprising, be- 
cause we expect a (negative) gravitational binding energy 
of the two black holes in the initial data, although there 
is also orbital kinetic energy and possible gravitational 
radiation in the initial slice not taken into account. 

We note that we have not used the traditional measure 



of mass for a black hole. 



MAH = yMi?, + ^ (5.1) 

which includes the angular momentum (and adds a small 
contribution to the mass), because the angular momen- 
tum of the final BH is not known in advance l46l| . 
Some previous studies 0,11 113 could safely assume that 
very little or no angular momentum was radiated, and 
were therefore able to use angular momentum of the 
spacetime to compute this mass, and then to calculate 
the total energy radiated. For the cases studied here, 
the angular momentum of the final black hole must be 
measured in each case, as we discuss below. 

Beyond simple determination of the AH masses, much 
more physics of the merger process can be obtained by 
studying the dynamics of the merged horizons. In this 
study, in all cases the black holes could be evolved to 
well past common apparent horizon formation, allowing 
numerous physical characteristics of the final black hole 
to be computed as it settles. 

We begin with a study of the shape of the horizons. 
First, we note that for evolutions that can be carried 
out to a point of near stationarity, when the final black 
hole has settled down, the event horizon can be accu- 
rately traced by integrating null surfaces backward in 
time through the spacetime [2^ i and further that a 
detailed analysis of the physics of these event horizons 
can be carried out [s^ . 14^ |4^ . We have carried out 
the analysis for these 3D spacetimes using the method 
described in In Fig. El we show a comparison of 

the coordinate locations of the event and apparent hori- 
zons for the near-ISCO case QC-1 at various times for 
a resolution of dx = 0.08M. On the initial slice t = 0, 
no common event or apparent horizon exists. Due to the 
small size of the horizons (compared to the grid spacing), 
the location of the event horizons can not be found with 
sufficient accuracy on the initial slice. For that reason we 
only show the location of the initial apparent horizons. 
By time t — 13M, the event horizon sections have just 
merged, with an interesting caustic structure that will be 
studied elsewhere. Interestingly, away from this caustic 
region the coordinate locations of the EH and AH agree 
with great precision, with the EH just outside the AH. 
By time t = 18M, we find the first appearance of a com- 
mon AH, which as expected occurs later than the merger 
of the EH. Finally, at the late time t = 30M, the AH and 
EH locations agree very well, indicating the final merged 
black hole has settled down. The close agreement in lo- 
cation of both local temporal measures (AH) and global 
temporal measures (EH) of horizon structures gives con- 
fidence in the overall accuracy of the merger evolution. 

The coordinate location of horizon structures are use- 
ful for comparisons and for topology measures, but they 
say nothing about pure geometric measurements from 
which physical quantities can be derived. As shown first 
by Smarr 49] for a Kerr BH, the spin parameter a/m 
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FIG. 4: Evolution of the irreducible mass Mirr of the AH as 
a function of time. Part (a) shows the mass for each QC-N 
model, for the low resolution dx = 0.08 evolutions, showing 
both individual- and common-horizon masses on the same 
scale. (For the QC-4 model, the individual AHs are briefly 
"lost" due to the proximity of the excision region, somewhat 
before a common horizon is found.) Part (b) shows the mass 
for the QC-0 model for the dx = 0.08M, 0.06M, and 0.048A'/ 
evolutions, with individual-horizon masses in the left sub-plot 
and common-horizon masses in the right sub-plot. For the 
dx — 0.08A/ QC-0 case, results from both apparent-horizon 
finders are shown, AHFinderDirect as a solid line, AHFinder 
as points. 



uniquely determines the geom etric shape of the BH hori- 
zon. As discussed in p5 l3Cll|. this relation allows one to 
use the horizon shape as a tool to measure the spin of a 
black hole in dynamic spacetimes (St |2j, [s^l • A reason- 
able measure of the shape of BHs that are nearly Kerr 
is the ratio of proper horizon circumferences measured 
along polar, Cp, and equatorial, Ce, geodesies (3^. The 
value of C'r — Cp/Ce is unity for a Schwarzschild BH, and 
less than unity for Kerr BHs, with the precise value de- 
pending uniquely on the rotation parameter a/m as dis- 
cussed below. As shown in 's^ , a distorted rotating BH 
should oscillate about its equilibrium Kerr shape. In ax- 




Event Horizons 
Apparent Horizons 

1 
X (M) 



FIG. 5: Event and apparent horizons for an evolution of the 
QC-1 model {dx = 0.08M), showing horizons on the initial 
slice {t = O.OM), when the EH first merges [t = 13. OM), at 
the appearance of the first common AH {t = 18. IM), and 
at a late time slice {t — 30. OM) when the system has essen- 
tially settled and the two contours cannot be distinguished 
anymore. 



isymmetry, the polar circumference is uniquely defined, 
but in full 3D one has multiple choices. In this study, 
we compute the ratio Cr along circumferii in the coordi- 
nate planes (because of symmetry about the plane of the 
orbit, the equatorial circumference is uniquely defined). 

Fig.Elshows results for the case just outside the ISCO, 
QC-1, plotting the evolution of the ratio of polar to equa- 
torial circumference Cr, measured in the xz plane, for 
both common apparent and event horizons, for the low 
resolution case with dx = 0.08M. The EH and AH 
shapes agree extremely well, as was also suggested by 
Fig. 1^1 Both curves reveal the oscillation pattern ex- 
pected in the geometry of a distorted rotating BH formed 
in the merger. 

Following 1^ [s^l , we use the oscillation pattern in Cr , 
and its equilibrium value, to extract important informa- 
tion. The final BH formed during a merger is expected to 
have its quasi-normal modes (QNMs) excited. For each 
member of the Kerr family, characterized by a/m, a set of 
discrete QNMs exists. Further, for each {£, m) pair, there 
is a fundamental mode and higher overtones. In princi- 
ple, all of these mode may be excited, but in practice the 
higher overtones are much more strongly damped than 
the fundamental. Although associated with the poten- 
tial barrier surrounding the BH, the QNMs excite grav- 
itational waves that propagate both to infinity and also 
inward, crossing the horizon where they manifest as os- 
cillations in the horizon geometry ^28^ _30J. Therefore, 
one should find various QNMs in the oscillations of Cr, 
and can use this information to determine the angular 
momentum of the final BH. 

The QNMs for rotating black holes have been studied 
extensively via perturbation theory, and have been com- 
puted as a function of a/m in Refs. [soUsij . Using these 
known values, and building on previous work l30j . we 
have devised the following procedure to analyze the hori- 
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FIG. 6: Comparison of EH to AH geometry for case QC-1. 
The ratio of the polar and equatorial horizon circumferences, 
Ct = Cp/Ce, after formation of common horizons, measured 
in the xz plane, is shown for the low resolution {dx = 0.08M) 
case. Both AH and EH horizons oscillate at the QNM fre- 
quency of the final Kerr BH, as confirmed by the fit shown 
(see text). 



zon oscillations. One can show analytically that for a 
Kerr black hole, the quantity Cr has a value given by 



^ _ 1 + y/l - (fl/m)2 ^ 



(a/m) 



(1 + - (aAn)2)V 

(5.2) 

where E{k) is the complete elliptic integral of the second 
kind 

/•V2 

E{k) = / VI - ksm^edO. 
Jo 

For simplicity Cr can be well approximated by the 
empirical-theoretical formula, accurate to within less 
than 1% Hi, 



- (a/m)2 + 1.55 



2.55 



(5.3) 



A perturbed BH horizon with a given value of a/m os- 
cillates about this equilibrium value, largely through a 
superposition of its QNMs [28L l30l |. Therefore, the func- 
tion Cr should be composed of an equilibrium value and 
a linear combination of various QNMs (G^l- We use an 
iterative least-squares fitting procedure to determine the 
value of a/m that best fits the extracted curve Cr • 

In these studies, we have found it is sufficient to con- 
sider only the offset, and the phase and amplitude of the 
fundamental modes (given by ^ — m) as fitting parame- 
ters. Here we have used the £ = 2,m = 0,1,2 modes, but 
we note that unless a/m becomes very large, the QNM 
frequencies are all very similar to each other. We begin 
with a fiducial set of modes for a given a/m value, adjust 
their phase and amplitude and the offset value, which de- 
termines the approximate value of a/m through Eq. (|5.3|) 



for the first iteration. We then use this value of a/m to 
choose the QNM frequencies for the next iteration, which 
through the offset determines the new value of a/m. For 
the data sets studied here, this procedure converges ex- 
tremely rapidly, independent of the starting value of a/m. 
We have also used a binary search technique over a/m 
to find a minimum error in the least squares procedure. 
This method gave identical final results for a/m. We note 
that the fitting procedures use the QNM frequencies for 
unit BH mass and vary only a/m. This is justified by the 
fact that Mbh ~ -^adm ~ 1- 

Turning back to Fig. for QC-1, we see that not 
only do the curves obtained for both the EH and AH 
closely agree, they both show the oscillations of the final 
Kerr BH formed in the coalescence process. The result 
of the fitting procedure to the Kerr modes, described 
above, is also displayed. We see a beautiful agreement 
between both the event and apparent horizon measure- 
ments, which are completely independent, and the fit to 
the QNMs of the final Kerr BH. 

This result shows conclusively that in the coalescence 
process, the QNMs of the final black hole are strongly 
excited, and dominate the evolution of the horizon after 
its formation. Although we have not attempted to ex- 
tract waves in the far zone in these studies, due to the 
nearness of the outer boundary, this result also implies 
that at late times QNMs will dominate the waveforms. 
Below, we study the determination of the spin of the final 
black hole from this procedure. 

In many cases we are able to evolve these systems sig- 
nificantly beyond the formation of the final BH, so that 
its properties and be studied in depth. We then have ad- 
ditional analysis tools at our disposal. For models where 
an event horizon was located (QC-0 and QC-1), the hori- 
zon generators can in principle provide an independent 
estimate of the spin of the final black hole, as shown 
in Isil Is^ . Several factors complicate the analysis 
of the horizon generators in these cases. First of all, for 
QC-0 and QC-1, the evolutions did last long enough past 
merger for us to track the EHs accurately during the 
merger, but not quite long enough to track them accu- 
rately at the time when the black hole has almost set- 
tled down to its final stationary configuration. It is only 
at this late time that the formula relating the generator 
angular velocity to a/m is correct. Additionally, the an- 
gular velocities are gauge dependent. Therefore we have 
to correct the coordinate angular velocities not only for 
the effects of our co-rotating coordinate system but also 
for any co-rotation shift developed by the F-driver shift 
(equations ^A.2\ and (|3.3|l ). We have, as yet, not found a 
satisfactory way of doing this and are consistently getting 
lower estimates for a/m using the generators compared 
to the results obtained by fitting QNMs to the EH cir- 
cumference ratios. 

The isolated horizon formalism [s^ Is^ l provides an ad- 
ditional tool to study the final horizon to extract phys- 
ical information. In this case we apply this formalism 
to the final merged AH to determine a spin by finding 
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dx Cr xz Cr yz Cr avg IH 

0.080 0.429 0.471 0.450 0.470 

0.060 0.548 0.597 0.573 0.570 

0.048 0.603 0.653 0.628 0.660 



TABLE II: Table comparing different angular momentum 
measures a/m for the final black hole in a QC-0 merger simu- 
lation at different resolutions dx. Cr refers to a/m determined 
via a QNM fit as described in the text. The cases xz, yz refer 
to the ratio of the circumference in the xz to the xy circumfer- 
ence and of the yz to the xy circumferences. The avg denotes 
the average of the two. IH denotes the value computed by the 
isolated horizon finder, measured on the plateau at t = 25A4". 



FIG. 7: Comparison of a/m measurements for case QC-0. We 
show the circumference ratio C,., measured on the AH, for 3 
different resolutions for the two different polar circumferences 
(computed in the xz and yz plane). We also show the results 
of the fitting procedure to determine an "offset" value of C- 
(straight lines) for the average of the measurements of circum- 
ference ratios as listed in Table HTl showing good agreement 
with the curves. The offsets determined from the two indi- 
vidual fits are used to compute the error bars at the left side 
of the figure. Also plotted is the result from the isolated hori- 
zon computation converted into a circumference measure for 
comparison. We note that the value from the isolated hori- 
zon has developed a plateau between t = 2Q — 25M, which 
agrees relatively well with the fits to Cr for the geometric 
measurements. 



an approximate solution to the Killing equation in the 
surface [s^. While such a solution might not exist im- 
mediately after the formation of the common AH we have 
found that within 7— lOAf the formalism provides reliable 
information as shown by the a/m measure approaching 
a constant value. 

We now apply both the horizon geometry and isolated 
horizon measurements discussed above to the case stud- 
ied in most depth for this paper, QC-0. First, as noted 
above, the geometric measurements Cr have been made 
in two different ways, using the polar circumference in 
both the xz and yz planes. In Fig. [3 we show the evo- 
lution of Cr in both of these planes, as well as a fit to 
determine the "offset" value, which allows us to estimate 
the spin of the final BH as described above. For com- 
parison, we have taken the measurement of the isolated 
horizon spin, and converted it into an instantaneous Cr 
via Eq. (|5.2|) . We see that this measurement settles into 
a plateau after about t = 20M, agreeing fairly well with 
the geometric measurements. 

In Table ITTl we turn these results into quantitative mea- 
surements of the spin of the final BH, for different meth- 
ods and different resolutions, all for the case QC-0. First, 
we note that at the lowest resolution, all measurements 
give consistently low values of the final spin, averaging 
to a/wiowres ~ 0.46. This estimate seems extremely low. 



dx 


0.080 


0.060 


0.048 


a/m 


0.450 ± 0.021 


0.572 ± 0.025 


0.632 ± 0.028 


Mirr 


0.947 


0.933 


0.923 


Mar 


0.973 ± 0.003 


0.978 ± 0.005 


0.980 ± 0.006 


Jrad (%) 


45.3 ± 2.9 


29.6 ±3.7 


22.1 ±4.5 


S.ad (%) 


3.61 ±0.25 


3.12 ±0.45 


2.97 ±0.59 



TABLE III: The table shows the radiated angular momentum 
and energy for the different resolutions dx for the case QC-0. 
Starting from Table ^ we use the minimum and maximum 
values for a/m at the different resolutions. The ranges given 
refer to these values for a/m. For the irreducible mass Mirr we 
take the values reported from the apparent horizon finder at 
t = 25M. We compute the horizon mass Mah which includes 
angular momentum and then compute radiated energy and 
angular momentum using initial data values of Madm = 1-01 
and (J/Af2)ADM = 0.779. 



however, for a system with initial J/M^ = 0.78. How- 
ever, the study of resolution effects on this quantity shows 
that it is very sensitive to resolution. At the highest res- 
olution studied, we measure the final BH to have spin 
a/m-hiros ~ 0.65. Unfortunately, the various measures 
used here do not show clear second order convergence 
properties, and therefore it is not possible to perform a 
Richardson extrapolation. 

Therefore, with present resolution and techniques, we 
can say with confidence that low resolution simulations 
underestimate the spin of the final black hole, and that 
at the highest resolutions possible in the present study, 
we believe we are still underestimating the final spin by 
a small amount. A separate investigation of 3D evolu- 
tion of single, spinning, distorted black holes is under- 
way to investigate resolution and computational param- 
eters needed to accurately resolve and determine the spin 
from horizon measurements |54l |. These results do indi- 
cate that the merger process, for the cases studied, ra- 
diates a small amount of angular momentum (less than 
about 20%). The results are summarized in Table Ulll 

Finally, we turn to the computation of the final black 
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hole mass. In Table Hill we also summarize computations 
of both the irreducible mass (not including spin) A/i^r 
and the horizon mass (including spin) Mah as defined 
in Eq. (|5.1|l 6^. In computing Mar we use a range of 
values for a/m, as listed in Table Hll to compute the error 
bars listed in Table ITTTl Note that there are additional 
errors in particular regarding the value of Mi„ which 
have not been included. We compute i^rad = (-^adm — 
Mah)/Madm and Jiad = [Jadm - ia/m)Ml^)/ Jadm- 
Here we see that, although we are unable to precisely 
determine the total energy radiated with present tech- 
niques, we nonetheless have demonstrated the high de- 
gree of accuracy that is now achievable for BH merger 
simulations. Our highest resolution simulations find a fi- 
nal BH mass within 3% of the ADM mass, and the trend 
as a function of resolution indicates the total energy ra- 
diated to be probably less than this. 

We note that Ref. '2&\ estimated, via completely differ- 
ent methods, that roughly 12% of the angular momen- 
tum is radiated, for the same initial data sets studied 
here, and that the total energy radiated is 3%. These re- 
sults are consistent with those that can now be obtained 
via full nonlinear evolution. These points are discussed 
further below. 



VI. DISCUSSION 

We have performed the first systematic study of binary 
black holes constructed to be in approximately quasi- 
circular orbits via effective potential methods, using the 
full machinery of 3D numerical relativity. The data sets 
chosen range from the innermost circular orbit, or ISCO 
(QC-0 in our notation), as determined from this construc- 
tion, to data sets with separations significantly beyond 
this, in what we refer to here as the near-ISCO regime. 
In each case, using a combination of recently developed 
techniques including excision, new gauges, the BSSN for- 
mulation, and a co-rotating frame, we have been able to 
carry out the evolutions until after a common apparent 
horizon formed, indicating an upper bound on the time 
of merger of the two black holes. 

The coordinate times at which the common hori- 
zon formed were found to be significantly less than the 
timescale of an orbit, defining this timescale from the ini- 
tial angular velocities. We also show that as one evolves 
initial data sets further out in this near ISCO regime, 
the holes merge after longer times but for progressively 
smaller fractions of an orbit. We compared these results 
with results from head-on collisions from corresponding 
distances, and found that the time to merger shows a 
very similar trend. This suggests that the initial data 
models studied here are already in an extended plunge, 
and not in quasi-circular orbits as one might expect from 
the effective potential considerations. We also found that 
the angular velocity of the co-rotating frame needed to 
keep the holes in a "head-on" collision as seen in coor- 
dinate space, is roughly 2/3 of what would be expected 



from considerations of the initial data, even after one 
takes into account the intrinsic co-rotating component of 
the shift condition used This result is consistent 

with our findings that the holes are merging much more 
rapidly than expected from initial data considerations. 

It is important to point out that these results are 
also consistent with two completely different studies, the 
Lazarus results about the near ISCO regime, and inves- 
tigations of binary neutron star systems near the ISCO. 

The same black hole data sets that we used here 
have been studied in th e hy brid perturbative-numerical 
Lazarus approach |3, 0, ■ This approach evolves the 
same data sets numerically, only to a "linearization" 
time, of roughly lOM, where the two holes are close 
enough that the system could be treated as a single 
perturbed BH, after which a perturbative approach was 
used to complete the evolution. Recall the results from 
Lazarus on the linearization time along the initial data 
sequence. According to Ref. ,25], for the cases studied 
here, QC-0 to QC-4, the linearization time increases with 
separation, but with a slope similar to, or even less than, 
that seen for coalescence times found here. This result 
gives support to the notion that this black hole sequence 
is not approaching an orbital configuration out to QC-4. 

Simulations of inspiral sequences of binary neutron 
stars were carefully studied by Miller and Miller et al. 
in [s^ Ist) ] , who examined the correlation of waveforms 
from model evolutions against PPN waveform templates. 
These authors find strong evidence that the circular orbit 
assumption can give rather poor approximations for two 
neutron stars in the near ISCO regime. Extrapolating 
these results to the binary BH case, note that the initial 
separations of QC-0 to QC-4 clearly fall into the regime 
where, based on [s^lSTj. one would not expect an orbital 
configuration. Again, this is consistent with our results 
that indicate a plunging configuration for all data sets 
studied here. 

A number of issues surrounding the merger time num- 
bers that we have presented need to be considered with 
care, in particular effects due to boundary treatment, nu- 
merical resolution, and the gauge choice. 

We have taken some care to assure that boundary ef- 
fects do not significantly alter the orbital dynamics and 
merger times picture. Boundaries have been placed at 
a distance causally separated from the common hori- 
zon, although even without this restriction we have ob- 
tained essentially identical runs with boundaries in both 
closer and further separations, indicating that they do 
not strongly influence the local dynamics of the black 
hole horizons for the timescales considered here. Though 
the numerical domain of dependence is larger, the Hamil- 
tonian constraint was found to converge to second order 
in the region of the common horizon. We also found that 
the merger times are not significantly affected by the res- 
olutions used in these studies. 

As a point of reference, in Ref. [s^l, high resolution 
(0.054 Af) was required to get convergent results. Our 
coarsest resolution of 0.08 Af is somewhat lower than 



11 



this, but our medium and high resolutions (0.06 M and 
0.048 M) are comparable. Unlike [s^], we were able to 
ensure that the outer boundary remains causally discon- 
nected from the merger region for all models studied. As 
a result, wc found that changes to the outer boundary lo- 
cation had only small effects on the time of merger (well 
within the error bars of Fig. |2J|. 

As remarked above, our studies indicate that for these 
BH simulations resolution and boundary location effects 
are not responsible for the quick coalescence times, espe- 
cially for the cases closer to the ISCO (e.g., QC-0,1). 

In addition to resolution and boundary effects, gauge 
effects will have an influence on the time at which a com- 
mon AH can first be found. In order to remove some 
potential slicing and coordinate ambiguity, for all cases 
we computed the proper time to common AH formation 
along the origin, which is at all times a point of symme- 
try between the holes. This measure of time is actually 
gauge invariant for EHs, but differs from time measured 
at infinity. The results show the same trend as above, 
that the holes are in an extended plunge and the greater 
separations do not show any tendency towards approach- 
ing an orbit. 

As described above, we have used a consistent set of 
gauge parameters for all the reported results (except the 
EH results). Variations of these parameters (to the ex- 
tent that was possible in order to still achieve a long- 
lasting run) were found to vary the time to appearance 
of a common AH by only a small amount, within the dis- 
played error bars. However, in the initial 5M of evolution 
when the black holes are still quite small on the grid, res- 
olution effects can influence the rate at which the gauges 
evolve (collapse of the lapse, for instance), resulting in 
an offset even with identically chosen gauge parameters. 
The error bars displayed in Fig. |21 are not statistically 
derived, but rather indicate the greatest range of results 
from many experiments with resolution and gauge pa- 
rameters. We have not, however, observed a signiflcant 
increase in the time to common apparent horizon forma- 
tion for any of the studied models. In particular for the 
QC-0 and QC-1 cases where a large number of runs under 
various conditions of gauge and resolution were carried 
out, we do not expect the merger time to vary greatly 
with changes to computational parameters (gauge, reso- 
lution, boundary). 

On the other hand, as indicated by the growing error 
bars, there is increasing uncertainty in the QC-3 and QC- 
4 cases. We note, for example, that the fisheye coordinate 
transformation that stretches the coordinates in physical 
space is very close to the initial puncture locations in 
these cases. This, or additional gauge effects, could con- 
ceivably affect orbital evolution properties. Further ex- 
periments for initial data in this regime and beyond are 
underway with mesh refinement, and with a wider range 
of gauges, and will be reported in future publications. 

We also applied a number of tools to measure various 
physical quantities from the final BH horizons formed in 
the merger, which were found to agree well. Both event 



and apparent horizons were determined for the QC-0 and 
QC-1 cases Geometric measures of the EH and AH 
horizon geometry were used to compute the mass and 
spin of the final BH, as well as the oscillation frequency 
of the radiation. In all cases, the final BH mass was 
found to be close (within 4%) to the ADM mass of the 
spacetime, providing a global measure of the accuracy of 
the simulations. Appearing with a delay of 4 — bM from 
the first joining of the EH, the common AH was found 
to coincide with the EH extremely well very soon after it 
merged. As we discuss below, in the highest resolution 
studies of the case QC-0, going beyond overall consis- 
tency of energy measures, we are nearly able to predict 
the total energy radiated. 

We also used the isolated horizon formalism to study 
the dynamics of the final BH, especially the final spin, 
and obtained similar results. It is interesting that within 
a few M of common AH formation, the estimate of the 
final angular momentum of the merged black hole from 
the AH distortion (via C^, defined above) and the same 
measure of EH distortion, agree well with the indepen- 
dent estimate obtained by the isolated horizon formalism. 
This suggests that it takes only a very short time for the 
final black hole to settle to a sufficiently Kerr-like state 
so that each of the approaches can be used. The fact that 
three very different measures of the spin of the final black 
hole agree so closely with each other gives us confidence 
in the vahdity of the results. 

However, the low resolution evolutions consistently 
under-predict the estimated spin by a significant amount 
for the cases studied. The resolution studies indicate that 
at higher resolution, the spin and mass measurements of 
the final BH horizon approach the asymptotic ADM val- 
ues, indicating that total radiation is small. For the QC-0 
case studied in depth, we estimate that roughly 20% or 
less of the total angular momentum is lost in the merger 
process. It is important to note that although the spin 
of the final black hole is sensitive to the resolution, the 
coalescence time is not. Even for the more distant config- 
urations in the pre-ISCO family studied here, resolution 
affects the coalescence time only relatively weakly, and 
the highest resolution studies find the holes to coalesce 
well before an orbital time period even for QC-4. 

We singled out the "ISCO" case, QC-0, for very de- 
tailed study. In addition to the calculations that indicate 
the system to radiate of order 20% or less of its angu- 
lar momentum, we also studied the radiated energy as a 
function of resolution. By taking into account the esti- 
mated angular momentum of the final black hole formed 
in the merger, we were able to compute the total mass of 
the final BH, Mah, defined by Eq. (|5.1|) . In all cases, it 
was found to be within 4% of the ADM mass, and at the 
highest resolution we found results consistent with 3% or 
less of total energy being radiated. 

Note that the results of the Lazarus approach 0, 
for the total energy and angular momentum radiated are 
consistent with our fully numerical results for the same 
initial data. The fact that the results are similar, yet the 
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evolutions were carried out with completely different evo- 
lution equations (ADM vs. BSSN) and gauge conditions 
(maximal/zero shift vs. l+log/Gamma-driver shift), and 
then were further differentiated by the linear and nonlin- 
ear approaches, gives added confidence in the results. 

The primary focus of this paper has been on the orbital 
dynamics of a family of binary black holes expected to be 
in quasi-circular orbits, and not on the waveforms gen- 
erated. However, by making detailed measurements of 
the apparent and event horizon geometries, we have been 
able to show conclusively that the quasi-normal modes of 
the final black hole are strongly and cleanly excited, and 
they are dominated by the fundamental I = 2 modes. 
These QNMs will therefore be very strongly seen in the 
waveforms far from the holes. 

While the simulations reported in this paper were be- 
ing analyzed for publication, another study was carried 
out using similar techniques, but with a different code 0. 
In that study, another binary black hole data set, from 
a family which is similar to the family studied here but 
with L/M ~ 9, further beyond the ISCO than the present 
study (QC-4 has L/M = 7.84), was found to evolve be- 
yond the predicted orbital time period (114M) without 
finding a common apparent horizon. Although our re- 
sults show a trend for this family of data sets to be 
in near plunge when L/M < 8, the results in Q for 
L/A'I — 9 are not necessarily in contradiction with the 
findings of this study. Further, there are a number of 
differences in the methodology, including gauge condi- 
tions and the computational grid setup. In this paper, 
we have stressed the importance of sufficient resolution 
to determine accurately the angular momentum of the 
final black hole. At the same time, we have not seen a 
strong dependence of the coalescence time on the res- 
olution for the range of simulations carried out here. 
We note that although our highest resolution runs are 
for dx — 0.04Af, slightly coarser than the typical fine 
grid resolution dx — 0.03125Af reported in we would 
not expect this alone to account for the differences. De- 
tailed study, beyond the scope of the present investiga- 
tion, will be required for data sets further into the near- 



ISCO regime to investigate these results. 

The question arises at which point of the QC sequence 
there is a transition from plunge to inspiral. For suf- 
ficiently large separations this type of initial data is 
known to be a good approximation to quasi-circular post- 
Newtonian data [Tolll3|. The results presented here in- 
dicate that the effective potential method does not pro- 
duce data sets that will evolve in quasi-circular orbits at 
initial proper separations less than 8M. We note that 
for a different initial data set the ISCO was found at 
an approximate distance of L/M = 6.61, much further 
out than with the method of minimum effective binding 
energy used for this study 

Recent progress in the use of mesh refinement in our 
codes 0, 113, as well as the introduction of higher or- 
der finite difference operators, will enable us to minimize 
resolution dependent effects. Work is in progress to inves- 
tigate the transition in the initial data between orbiting 
and plunging black holes, and to study the sequence of 
initial data sets developed usin g th e HKV conformal thin 
sandwich approach p!4[l5l[l7ll59|. 
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